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Abstract 

We show that non-perturbative Gxed points of the exact renormalization group, their 
perturbations and corresponding massive Geld theories can all be determined directly in 
the continuum - without using bare actions or any tuning procedure. As an example, we 
estimate the universal couplings of the non-perturbative three-dimensional one-component 
massive scalar field theory in the Ising model universality class, by using a derivative 
expansion (and no other approximation). These are compared to the recent results from 
other methods. At order derivative-squared approximation, the four-point coupling at zero 
momentum is better determined by other methods, but factoring this out appropriately, all 
our other results are in very close agreement with the most powerful of these methods. In 
addition we provide for the first time, estimates of the n-point couplings at zero momentum, 
with n = 12, 14, and the order momentum-squared parts with n = 2 • • • 10. 
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1. Introduction 



The motivations for the present paper are three-fold. Firstly, we wish to expand upon 
the renormalization group reasons, given in ref.[Q, for the quantization of renormalised 
couplings. Since these couplings correspond to relevant (and marginally relevant) pertur- 
bations about a fixed point, this necessitates first briefly reviewing how the renormalization 
group also 'self-determines' the fixed point structure itself -[Q: as we recently empha- 
sised 0, the fixed-point and its eigen-operator spectrum - and hence the corresponding 
massless quantum field theory, can be deduced and computed entirely from consistency 
arguments applied to the effective action and its exact renormalization group flow close 
to the fixed point - i.e. without needing to go through the construction of introducing 
an overall cutoff Ao, a sufficiently general bare action, and then taking the continuum 
limit Aq — > 0. Thus the universal continuum properties are accessed directly without this 
standard, but for quantum field theory, actually artificial and extraneous, scaffolding. 

It would seem rather strange if having been able in this way to derive directly the 
continuum massless theory, the corresponding massive theory required the reintroduction of 
the concepts of bare couplings and corresponding tunings to reach the continuum limit. At 
first sight this appears to be the case, however the second and main purpose of the present 
paper is to show that once the massless theory's fixed point and eigen-operator spectrum 
are known, the massive theory may be constructed, again directly, without recourse to 
any limiting procedure. An important byproduct of this analysis, is the transformation 
of the continuum limit of the Wilsonian RG (Renormalization Group) into the form of a 
self-similar flow for the underlying relevant couplings, in close analogy to the usual field 
theory perturbative RG, although here the (3 functions are defined non-perturbatively.0 

We concentrate on a description within the derivative expansion approximation to the 
renormalization group[|J (see also fl9fl-||ll|), but point out at the appropriate points how 
and why precisely the same effects can be expected to work in determining the fixed points, 
eigen-operator spectrum and massive continuum limits also for the exact renormalization 
group. 

Finally, we apply these concepts to a calculation, which however is interesting in its own 
right: the universal coupling constant ratios of the three dimensional Ising universality class 
scaling equation of state. In quantum field theory terms, this corresponds to determining 
the one-particle irreducible n-point functions at zero momentum (equivalently the O(p ) 



1 Some related remarks can be found in refs.||— fl. 
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terms in a series expansion in powers of the momenta) for the three-dimensional one- 
component massive scalar field theory based about the non-perturbative Wilson-Fisher 
fixed point. (These are universal once they are divided by appropriate powers of the scalar's 
mass.) We compute these universal ratios both for the Local Potential Approximation of 



the Wegner- Houghton equation |T!| |13|], and for a smooth cutoff as utilised in refs.[||||. 
The latter allows also to go to next order - 0(d 2 ) - in the derivative expansion. We will 
see that the results improve, and compare well with the very recent high order perturbation 
theory results of Guida and Zinn- Justin [|14| . Indeed, appropriate universal higher coupling 
constant ratios lie within the errors of the corresponding resummed perturbation theory 
results|n[]. Moreover, we are able to give estimates for even higher point couplings (n = 
12 •••14), and also to estimate the corresponding universal 0(p 2 ) terms in the first six 
n-point functions. At present, these are not available from any other method. 

Computing the n-point couplings of the scaling equation of state in this way is actually 
a slightly eccentric use of the derivative expansion approximation, since the results most 
naturally present themselves as full numerical solutions for the equation of state. Indeed 
in ref. [|15|, a numerical computation of the scaling equation of stateH has already been 
presented. These authors obtain the equation of state by tuning an appropriate bare ac- 
tion, and solving the RG within a modified O(d ) approximation incorporating anomalous 
scaling (viz. r/). Our purpose for computing instead the corresponding Taylor expansion 
coefficients, was primarily to provide a direct comparison of the first two orders of the 
derivative expansion per se with the emerging accurate resummed perturbation theory re- 
sults, preliminarily announced by Zinn- Justin, and some independently also by Sokolov ||17|| 
at the August conference "RG96" in Dubna. For this comparison, it is sufficient to com- 
pute only in the symmetric phase. (The conceptual points of the paper however apply 
to either phase.) A similar computation in the broken phase would be very interesting 
and allow estimation also of 'amplitude ratios' |HJ by the derivative expansion. Alterna- 
tively these could be computed from the couplings obtained here, by analytic continuation 



in the mass ||14|| . Extensions to three dimensional O(N) symmetric iV-component scalar 
field theory, and to massive two dimensional theories, in particular the infinite sequence of 
multicritical fixed points, described by one-component scalar field theory, seems straight- 
forward, the formalism to 0(d 2 ) for these cases having already been developed HP|[|l8|. 



Also, it is certainly possible to estimate even higher point couplings, using a more serious 



See ref. [16] for a lattice field theory computation. 
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numerical attack than we attempt in this paper. (All the computations for this paper were 
performed within the Maple package.) 

A more detailed precis of the present paper now follows. In sec. 2, we set up some of the 
notation and recall some salient facts: the extraction of the Legendre effective potential 
(and action) from the A — > limit of the Wilson effective action - which is a necessary 
step to construction of the equation of state from the exact RG, the self-determination 
of the fixed point structure through the requirement of non-singularity, the large field 
behaviour and thus determination of the fixed point's Legendre effective potential. We 
point out that the leading large field asymptotics are independent of the effective cutoff A, 
physically a reflection of the fact that quantum corrections to the free energy are negligible 
for sufficiently large fields. 

In sec. 3, we study the spectrum of perturbations about the fixed point (a.k.a. eigen- 
operators). We show that quantization of the spectrum of dimensions of these operators 
follows from the requirement that the scale dependence can be absorbed in an associated 
coupling, i.e. a renormalised coupling, and thus correspond to the universal self-similar 
flow characteristic of the continuum limit. This is not guaranteed beyond the linearised 
level. By examining the flow exactly, which can be done for large fields, we establish that 
only perturbations behaving as a power of <f> for large field (f>, can be associated with a 
renormalised coupling. As a byproduct we derive the form of the Legendre effective poten- 
tial to first order in the associated physical coupling. While these results are established 
within the Local Potential Approximation (LPA), we explain why (as in sec.l) the same re- 
sults should follow also for the exact RG. Finally in this section, we apply Sturm-Liouville 
theory to the LPA to prove that those perturbations that do not behave as power-law for 
large <fi, collapse under evolution with scale into an infinite sum of perturbations that do 
behave as power-law. Thus all continuum physics is described in terms of the fixed point, 
the 'power-law perturbations', and their associated renormalised couplings. The main 
statements of this section confirm and generalise the corresponding findings of ref. |l| . 

The previous two sections thus describe how to obtain directly any massless continuum 
limit, and its associated operator spectrum, "directly" because no limit is actually taken. 
In sec. 4, we determine how to set up directly the corresponding massive continuum limits. 
We first review the standard lore on how a non-perturbative massive continuum limit 
is obtained, defining in the process the so-called renormalised trajectory. We show that 
this trajectory can be defined directly in terms of a boundary condition on the flow, 
involving a coupling g, such that finite values of g yield finite continuum limits, and 
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obtain as a consequence the leading terms in large <fi dependence of the corresponding 
massive continuum limit. The boundary condition still involves a limiting procedure (its 
dependence on A is prescribed as A — > oo), but we show that this limit can also be 
removed by reexpressing the flow in terms of the underlying renormalised couplings and 
their associated beta-functions. In addition we set up a change of variables so that the 
universal ratios of couplings (viz. physical scale independent ratios) are obtained once all 
modes have been integrated out (i.e. as A — ► 0). Again in this section, we point out how 
and why these observations should hold for the exact RG. 

In sec. 5, we formulate the theory of sec. 4 specifically in terms of the LPA of the 
sharp cutoff (a.k.a. Wegner-Houghton) flow equation, applied to three dimensional one- 
component scalar field theory governed by the non-perturbative Wilson-Fisher fixed point, 
particularly the universal coupling constant ratios of the Legendre effective potential 
(equivalently scaling equation of state). We discuss the numerical implementation of the 
boundary conditions and flows, expansion directly in terms of the (ratios of) couplings, 
and the corresponding choice of (and independence of) closure ansatz. The short sec. 6 
delineates the differences that arise in using the 0(d ) approximation, equivalently LPA, 
of the smooth cutoff flow equation derived in ref. . 

In sec. 7 we go beyond the LPA, to order derivative-squared, thereby also allowing esti- 
mation of universal ratios of the momentum-squared terms in n-point Green functions, and 
more accurate estimation of the scaling equation of state. Wave-function renormalisation 
and running anomalous dimensions are dealt with explicitly here. Once again we provide 
here also the explicit formulae needed for the numerics. 

Finally, in sec. 8 we present the numerical results. We present numerical evidence show- 
ing independence of closure ansatz, for three different choices of ansatz. Picking the most 
rapidly converging ansatz, we deduce the sharp cutoff LPA, the O(d ), and 0(d 2 ) results 
for the universal coupling constant ratios, and combine the derivative expansion results 
into one final number together with an estimated error of truncation of the derivative 
expansion. These are compared to results from resummed perturbation theory [fT4"| . The 
perturbative methods are more powerful for low order couplings but the derivative ex- 
pansion eventually wins out for higher order couplings. Following ref.jnj], we then factor 
out an overall normalization (effectively the size of the four-point coupling) and present 
derivative expansion results (and sharp cutoff LPA results) for the corresponding ratios 
F21-1. These are seen to be much more accurately determined by the derivative expansion, 
and are in close agreement with the most accurate determinations from other methods. 
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(A table of comparisons is presented, including the results from resummed perturbation 



theory, e expansion, the exact RG approximation of ref.|15|, high temperature series and 
Monte-Carlo estimates.) We thus conclude that the main error in the derivative expan- 
sion is in determining the lowest point coupling(s), and can be absorbed in an effective 
normalisation factor Ceff- 

There is no need for a separate summary and conclusions, since this has already in 
effect been incorporated in the above introduction and precis. 



2. Fixed Points 



First, let us review how the fixed points are determined by consistency arguments 
alone. To do this, we find it easier to concentrate on a specific equation and then to 
point out how the arguments generalise. Thus consider the so-called LPA (Local Potential 
Approximation) |13| to the Wegner- Houghton renormalization group [|12|. We remind the 
reader that the latter is just0 Q the sharp cutofii limit of Polchinski's form[[L9| of Wilson's 
continuous RGpOfl (Renormalization Group): 



dS A 

dA 



1 dA uv f 5S A 5S A 5 2 S A , x 
2 tr dA 1 56 56 ^ 1 ^ 



SSa ] 
56 J 



(where Auv(q, A) = Cuv/q 2 is a cutoff massless propagator. For more details on this 
equation see refs. || |L9| J?(J . It will not however be used in the rest of the paper.) The 
approximation corresponds to projecting on an effective action of the form 



S A = Jd D X {^(9 M </>phys) 2 + Vphys(</>phys, A)} 



i.e. to discarding all momentum dependent terms generated by the RG. (D is the space- 



time dimension.) For a single component scalar field, the result is ||13||||21||||22|| [Plp I 



JV(0, t) + d6V'(6, t) - DV(6, t) = In [1 + V"(6, t)] 



(2.2) 



Here t = ln(// p /A), \x v is an arbitrary physical mass scale, the primes stand for differentials 
with respect to 6, and d is the scaling dimension of the field 6. Generally this is given at a 



3 We look first at this long established version of the LPA. Later we will employ smooth cutoff 
versions, which also turn out to be more accurate. The subtleties involved in handling a sharp 



cutoff, both technical and physical, have been extensively addressed in ref. [11]. 
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fixed point by d = \ (D — 2 + rj) , however a consequence of the LPA is that the anomalous 
dimension r\ = 0. In this equation the scaling dimensions of the field and potential have 
been used to construct dimensionless equivalents Q: 

<j> = <p P /A d and V(</>,t) = V p (<j> p ,A)/A D . (2.3) 

Thus p-subscripted quantities refer to dimensionful physical quantities. This step is equiv- 
alent to (but more direct than) the traditional rescaling of the cutoff back to its original 
size after a blocking[p0|p^1. 

Actually, a further purely numerical transformation is necessary to obtain the properly 
normalised physical quantities appearing in (|2.1|): 

4> P = C^phys and V p = ( 2 V phys , (2.4) 

where C = (^) D / 4 y/T(D/2) was chosen to prettify the eqn. ( |2.2|) . This step will only need 
to be taken into account in the final sections containing the numerics. 

It is worth recalling here that the effective potential of the Wilsonian effective action 
S\ tends, in the limit A — > 0, to the Legendre effective potential ||@ (i.e. generator of one- 
particle irreducible Green functions evaluated at zero momentum). The latter is of course 
precisely what we will need later, to construct the equation of state. This is one conse- 
quence of the deeper connection also proved in ref.0, namely that the Legendre effective 
action (minus its infrared regulated kinetic term|J]): F\, for a field theory with infrared 
cutoff Cm(q,A) imposed, is a generalised Legendre transform of Polchinski's Wilsonian 
effective action S\ computed with ultra-violet cutoff Cuv^Q-, A) = 1 — Cm(q, A): I 

S A [<j ) } = r A [<p c } + ±(<j ) c - ( f ) ).-^.(<p c -<p) . (2.5) 

<^IR 



(Somewhat similar remarks have been made by other authors ]24J] ||25|| . N.B. </> c in this 
equation is the "classical field" , defined as a function of (p through the Legendre transform, 
in the usual way, i.e. by differentiating ( |2.5| ) by <fi at fixed </> c .) In this way we can have 
our cake and eat it: fixed points in the Wilsonian RG have a natural physical explanation 
through thinking about blocking, but at the same time we can recover the information 
apparently blocked away - via this generalised Legendre transform relation. 



4 providing that -^Cuv\q=o = 0, as follows trivially from observations in refs.[p|] pi] [j23[| 

5 It does not matter here whether quantities are physical, or dimensionless combinations. 
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At a fixed point, V(<j>, t) = V(<j>) is fixed i.e. independent of t. From eqn. (pT2|) , it 
would appear at first glance that there is a two parameter set of fixed point solutions V((j)) 
(because (|2.2| ) reduces to a second order ordinary differential equation), but this is true 
only locally. In fact all but typically a discrete number! of these solutions are singular 
at some finite value of the field <f> @ || [0 1| . Clearly the presence of such a singularity is 
unacceptable e.g. in leading to violations of Griffiths' analyticity. From our experience, the 
solutions that are well-defined for all real (p (typically then countable in number) , may all 
be identified with approximations to the expected fixed points of the exact RG, in the sense 
that they have the right qualitative properties and yield reasonable quantitative answers 
|| || [U • Importantly there appear to be none of the 'spurious fixed points' that afflict 
higher order truncations in powers of the field 

The form of the fixed point solution for large <p may be ascertained from eqn. (|2.2|) H : 

V(<j>) ~ A<j) 1+S as oo , (2.6) 

with A > an as yet unkown coefficient. We have used a hyperscaling relation to write 
D/d — 1 as the critical exponent 6. (We are assuming d > here and from now on. 
This thus excludes the two-dimensional exceptions 0. All scaling relations trivially hold 
automatically in derivative expansion approximations since scale invariance is not broken 
and the massless limit may be directly considered, as we are indeed so doing here.) The 
leading term above remains valid beyond the LPA and for any cutoff in the form wri tteS 
because it simply arises from solving the left hand side of eq. (|2.2|) ,i which in turn is purely 



generated by the assignment of scaling dimensions d and D to the field and potential 
respectively. This power law is thus precisely as required so that V p (<p p ,A) has a finite 
non- vanishing limit as A — > 0. Indeed, from ( |2.3| ) and ( |2.6|) , we obtain 



V p (c/> p ,0) = AcPl+ 5 . (2.7) 

In view of the relation to the Legendre effective potential outlined above, this is nothing 
but the free energy per unit volume, at the critical point, as a function of the field ||. 
Notice that the behaviour of V p ((j) p ,A) in the regime <\> p » A d , is also given by fl2.7j) , 
and in particular is thus actually independent of A. Physically, this is precisely to be 



6 an exception being critical sine-Gordon models in two dimensions || 

7 i.e. with 5 = D/d — 1, and d the full scaling dimension following from a non-zero r\. 

8 the neglect of the right hand side being justified by inspection 
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expected: it reflects the fact that quantum corrections to the free energy are negligible for 
sufficiently large fields. Indeed the asymptotic behaviour ( |2.6| ) for a fixed point, is the only 
such behaviour consistent with this 'mean-field-like' evolution i.e. one in which quantum 
corrections are absent. 

It seems reasonable to suppose that the same characteristics, as outlined above, also 
completely determine the fixed points in the exact renormalisation group. In this case@, 
the fixed point equation is a second order non-linear functional differential equation and 
presumably, locally has a full functional space worth of solutions but globally only typically 
a discrete number which are well defined for all 0(x), these latter being again the only 
bona fide fixed points. 



3. The eigen-operators 

To obtain the spectrum of (integrated scalar) eigen-operators we linearize around the 
fixed point by writing Sa[</>] = S[(f>] + ee xt s [<j>], or in terms of the LPA: 

V(<j )l t) = V( ( p)+ee xt v(<P) , (3-1) 

to first order in e. Here we have used the fact that the RG is quasilinear in t, and used 
separation of variables. Thus we have from (|2 . 2| ) : 

Xuffl + - DvM = 1 • ( 3 - 2 ) 

This time the ordinary second order differential equation is linear with non-singular co- 
efficients and therefore we are guaranteed to have for every real value of A, a continuous 
one-parameter (up to the arbitrary normalisation) set of globally well-defined real solu- 
tions. How can we square this with the fact that experiment, simulation etc. , typically 
only uncover a discrete spectrum of such operators and in particular only a discrete num- 
ber of relevant directions corresponding to eigenvalues A > 0? The answer was given in 
ref. |]] || (if rather compactly): only the discrete set of solutions for A and v ((/)), where v((f>) 
behaves as a power of cf> for large field, can be associated with a corresponding renormalised 
coupling g(t) and thus universal self-similar flow close to the fixed point. 

This comes about as follows. First note that if v((p) is to behave as power-law for large 
<fi, then from ( |3.2| ), this power is determined as 

v(4>) ~ acj) 1+5 - x/d as <j> oo , (3.3) 
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with a an as yet unkown coefficient. Again, as with (|2.6| ), this remains valid beyond the LPA 
and for any cutoff since the power is determined only by the left hand side of (|3l|) . Once 
this power law is imposed for </> — > oo, and for <fi — > — oo with possibly different coefficient,i 
and linearity of ( |3.2|) taken into account, we see that we have sufficient boundary conditions 
to overconstrain ( |3.2| ) leading to typically a discrete set of solutions for the eigen-values A 
and eigen-operators v((f>) 0. Therefore we need now to explain only why we require v(4>) 
to behave as a power-law for large <p. 

Studying eqn.( |3.2| ), and using (|2.6|), we see that if v ((f)) does not behave as ( |3.3| ) then 
for A > 0, it must behave as 

v(</>) ~ expA(D - d)(p 1+s as -+ oo , (3.4) 

where we are neglecting a multiplicative power-law correction. Only for the Gaussian 
fixed point V(<f>) = 0, does A = 0, and in this case we have instead v((f)) ~ expd(f> 2 /2. The 
behaviour (|3.4| ) is changed by going beyond the LPA, and depends on the choice of cutoff0 
because it results from balancing the right hand side of (|3.2[) against the derivative term 
of the left hand side. The crucial point for us however is that if v((f>) does not behave as 
( P-3|) , then it diverges faster than power-law as (p ~^ oo. And this statement is presumably 
universally true, holding for any (valid) choice of cutoff and also for the exact RG. 

From eqn. (|3.1| ), it is tempting to identify g(t) = ee xt as a (renormalised) coupling 
conjugate to the operator v((f>), but this t evolution followed from linearising in e, and 
for any finite coupling, no matter how small, there is a regime of large fields where such 
linearisation is not justified. [Thus e.g. ev(4>) » 1 for » (— lne) d ^ D in case (|3.4|) .1 



In this regime however, it is straightforward to solve exactly for the t evolution because 
the right hand side of (|2.2|) contributes negligably compared to the at-least-power-law 
behaviour of the left hand side, while the left hand side -being linear- is easily solved. 
Thus starting at t = with V(<j>, 0) = V((f>) + ev((f)), we have 

V((f),t) ~ e Dt V((f)e- dt ,0) =V((f))+ee Dt v((f)e- dt ), as </> -> oo , (3.5) 



where we have used ( |2.6|) (and are assuming, here and later, that subleading corrections 
to each term can be neglected). Indeed, once again, the neglect of the right hand side of 



9 or a boundary condition following from symmetry considerations, is imposed at the origin| 
10 e.g. compare the formulae of ref.|l. 
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( p.2|) is just the statement that quantum corrections are negligible in the regime of large 
fields; combining ( |3.5| ) and (|2~3j), 

V p ((f) p ,A)~^V(<t) p /^,0) for </> p »A d . (3.6) 

We expect then that ( |3.5| ) and ( |3.6| ) hold also for the exact RG. Now, for the power-law 
perturbations ( |3.3|) , (|3.5|) simplifies to 

V(</>, t) ~ A(P 1+S +ee xt a<J) 1+s - x / d as -»• oo , (3.7) 

which is just the large ^> behaviour of ( |3.1|) - confirming the identification g(t) = ee xt . 
Correspondingly from ( |3.6|) : 

V p (<f> p , 0) = A</> p +<5 + en* a 4> P +5 - X/d , (3.8) 

from which we see that g{t) corresponds to a dimensionful /i p -dependent physical coupling 
g p = e/ip of dimension A. But for the non-power-law perturbations ( |3.4| ), and in general for 
perturbations rising faster than power-law, the t dependence in (|3.5|) cannot be combined 
with e and thus absorbed into a running coupling. (Nor can all the \i p dependence be 
absorbed in a physical coupling.) We have thus confirmed the italicized statement at the 
beginning of this section, and expect that the statement holds beyond the LPA and indeed 
also for the exact RG. 

Finally, we know from Sturm-Liouville theory [p6f applied to (|3.2| ) that as soon as t > 0, 
non-power-law perturbations [i.e. behaving as ( |3.4j) , fl3~5|)1, can be reexpanded in terms of 
a sum over the power-law perturbations. It follows then that in any case only these are 
needed to span all the continuum physics. 

To prove these statements, we note that the differential operator "D, defined by writing 
( jOl) as Vv{4>) = A v ((/>), is Hermitian with respect to the weight function 

p(<P) = {l + V"(<f>)}exp-d y#^{l + y"(^)} , (3.9) 

and the power-law boundary conditions ( |3T3} ) (imposed at <p = ±oo). This follows because 
the large <fi behaviour of p is sufficient to ensure that the bilinear concomitant (the boundary 
term generated by integration by parts) vanishes. Indeed from ( |2.6|) , 

p{(j>) ~ exp-A(D -d)(j) 1+s as -»• oo , (3.10) 

11 We checked that at subleading order in the asymptotic expansion, the t dependence can still 
be entirely absorbed in g(t). 
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for A > 0, otherwise A = and p ~ exp — d(p 2 /2 (neglecting again, here and in the 
ensuing, the unimportant multiplicative power-law corrections). Now, from the general 
Sturm-Liouville theory, we have that the power law eigen- functions Vi((f>) form a discrete 
set, i = 1, 2, ■ ■ -, with corresponding, possibly finitely degenerate, eigen-values Aj, such 
that there is a most positive eigenvalue, which we call Ai, and an infinite tower of negative 
eigenvalues (i.e. taking the Aj's to be non-increasing in i, we have Ai — ► — oo as i — ► oo. 
Note that these facts concur with physics expectations) . The eigen- functions may be taken 
to be ortho-normalised: J_ d(f)p((f)) Vi((f))vj(4>) = From (|3.4f) and below it, and ( |3.5|) , 
we have that a non-power-law perturbation v((f>,t) = V((f),t) — V(<p) behaves as 

A>0: v(<f>,t) ~expA(D-d)e- Dt <j) 1+s 

as 4> — > oo , 

A = 0: ~ expde- 2dt (f) 2 /2 



thus the expansion coefficients Ci(t) = J_ d<pp((f)) v (</>, t)vi(4>), are well defined for all 
^ (A > 0, or t > ±§| 



t > 0, while for t > Hj- (A > 0, or t > ^ with A = 0), we have the completeness relation 



J d(t>p{<j>) lv(<f>,t) -Y^Ci(t)vi((/)) j> ^0 as iV 



■OC' 



We have thus proved the italicized statements of the previous paragraph. Note that these 
statements are the generalisation to any fixed point, of the corresponding statements made 
about Gaussian fixed points in ref. M . 



4. Massive continuum limits 

From this detailed study we actually have in place everything we need to set up directly 
the massive continuum limit, i.e. without actually having to perform the limit. Before 
doing this however, let us recall the standard lore|2(| on how a non-perturbative massive 



continuum limit is obtained. This is illustrated in fig.l. In the infinite dimensional space 
of bare actions, there is the so-called critical manifold, which consists of all bare actions 
yielding a given massless continuum limit. Any point on this manifold - i.e. any such 
bare action - flows under a given RG towards its fixed point; local to the fixed point, the 
critical manifold is spanned by the infinite set of irrelevant operators. Let us assume that 
emanating out of the fixed point there is just one relevant direction; the generalisation 
to more than one is straightforward and will be indicated later. Choosing an appropriate 
parametrization of the bare action, we move a little bit away from the critical manifold. 
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Fig.l. Tuning to a massive continuum limit, x marks points on the critical manifold, 
whereas the black blob is a point shifted slightly off the critical manifold. 

The trajectory of the RG will to begin with, move towards the fixed point, but then 
shoot away in the relevant direction (exponentially fast in t) towards the so-called high 
temperature fixed point which represents an infinitely massive quantum field theory. 

To obtain the continuum limit, and thus finite masses, one must now tune the bare 
action back towards the critical manifold and at the same time, reexpress late t physical 
quantities in renormalised terms appropriate for the diverging correlation length. In the 
limit that the bare action touches the critical manifold, the RG trajectory splits into two: 
a part that goes right into the fixed point, and a second part that emanates from the 
fixed point out along the relevant direction. This path is known as the Renormalised 
Trajectory [2(J (RT); the effective actions on this path are 'perfect actions' |27j] . In terms of 
renormalised quantities, the large t section of this path obtains a finite limit, namely the 
effective action of the continuum quantum field theory. 

Therefore, to obtain the limits directly we must first describe the RT. This is given, 
close to the fixed point, by 

V(<t>,t)=V(<t>)+g(t)v(<t>) , (4.1) 
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where the relevant coupling g(t), corresponding to the one relevant direction, satisfies 

g(t) = ge xt , (4.2) 

with g a constant. Eqn. ( fO| ) is the unique solution that 'tracks back' into the fixed point 
as t — > — oo. It satisfies the flow equation there [by linearisation as in ( ft.l| )1 since g(t) — > 0. 
Using (|4.1| ) as a t = — oo 'boundary condition' for the flow equation, we may now describe 
the whole RT. 

In a case where there is more than one relevant direction, we then of course have a 
continuum of possible RT's, which may be labelled by the ratios of the relevant couplings, 
thus labelling different physics for given overall mass-scale (a.k.a. correlation length). 
Eqn. ( f4Tl~D is replaced by a sum over the relevant couplings, each of which satisfies an 
equation of the form (|4.2) as t — > — oo. These generalizations are straightforward and thus 



will not be considered further. 

Note now, that any finite value of g yields a finite limit for the continuum quantum 
field theory. To see this, note that since A > 0, comparison of ( |3.3| ) with ( [2.6| ) shows that 
for any finite t, there is a sufficiently large <p where the linearised solution ( }4.1| ) is still valid! 
Indeed, the linearisation (|3.2|) is valid for (|4.1|) , providing we have g(t)v" ((/))/ V" (</>) « 1, 
and this is true for all <f> » Cg^e dt . (Here C = °^ 1+ ^a/i+($)~ A ^ > an< ^ we nave used the 
hyper-scaling relation (3 = dv = d/X.) Taking the limits t —>■ oo and — > oo, while obeying 
the above inequality, we may continue to use ( |4.1| ) and ( |4.2| ) and thus, similarly to (|3.7| ) 
and (|3.8|) , we obtain the scaling free energy per unit volume for large constant fields, for 
the system perturbed from the fixed point: 

W P ,O)~0j+* {A + ag^-W} for 4> p »Cg? . (4.3) 

Here we have written g p = gfip. The limit is finite for any finite g, as claimed, at least 
in this large <p p region, and we see that it corresponds to a finite physical coupling g p . 
Furthermore, if V p (<p p , 0) is finite for all 4> p , we have using (|2.3|), that 

V(M~ Vp (AU,0)/A D (4.4) 

should satisfy ( |2.2| ) in the limit A — > 0. Substituting, we see that ( |2.2| ) is indeed satisfied, 
in fact up to corrections ~ ln(A)/A D . 

We remark that higher order corrections to ( f4.3|) are presumably directly calculable by 
continuing the perturbative expansion fl4.1| ) to higher order in g, this being determined 
uniquely by similar considerations to sect. 3. 
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Obviously the form of the boundary condition: (|4.1|) , with (O) and t — > — oo, is rather 



inconvenient, indeed is yet another limit, and the flow equation is numerically unstable 
(stiff) in such a region - since perturbations behave exponentially in t. We can resolve 
both these problems by expressing the continuum RG in terms of self-similar flow of the 
renormalised couplings (i.e. analogously to the field theory perturbative RG). In other 
words we rewrite in this case, ( |2.2|) as a differential equation for V (</>, g) together with 
a (non-perturbative) beta function (3(g) for g(t). The advantage of this is that the flow 
equation in g has perturbations that behave only as a power-law (around the fixed point) 
and, as we will see, has remarkably good numerical behaviour. At the same time we have 
resolved the limit since the boundary condition defining the RT, may now simply be stated 
as0 

■?-V(<t>,g) =v{cj>) . 
dg g =o 

Rather than furnish further details at this stage, we combine this idea with another change 
of variables. Since there is really only one dimensionful parameter g p ^ it merely serves 
to fix the physical mass scale m p . If we provide a physical definition for this mass-scale 
(equivalent to a renormalization scheme for g p ) then the dimensionless coupling constant 
ratios formed using powers of m p , will be universal. Evidently these are generated by the 
dimensionless potential (and field): 

U(<p) = V p {m}<p,0)/m° , (4.5) 

from which it follows that the equation 

U'(<p)=h , 

is nothing but the universal scaling equation of state (h being the external 'magnetic field', 
equivalent to the source J in quantum field theory). Let us define a (dimensionless) running 
'mass' m = m(g) with the property that it corresponds in the limit to the physical mass 
m p , i.e. 

m p = hm Am(g) . (4-6) 
It will be helpful to change variables from g to m. Indeed, defining 

U(<p,m) = V{m d <p,t)/m D , (4.7) 



12 Similar remarks have been made in ref. ftl. 

13 corresponding to the one relevant direction, all fi p dependence being absorbed in g p 
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and using (|2.3| ) , ( fOj) and ( 4.5|) , we see that U (ip, m) has the property that it tends to 
the universal effective potential ( f4.5[ ) in the limit m — > oo. By rewriting the RG ( |2.2j ) 
as a flow equation for U((p,m), we obtain the universal physics directly as the large to 
behaviour. The flow with respect to t is never needed explicitly, having been 'hidden' in 
the beta-function for m: dm/dt = /3(m). 

In the next three sections we flesh out the above sketch of changes of variables for the 
cases of LPA with sharp cutoff, LPA - equivalently O(d ) - for the smooth cutoff employed 
in H, and 0(d 2 ) for this smooth cutoff. In all cases we study only the symmetric phase 
and use as definition of m p , 

7^ = ^(0,0) , (4.8) 
in concord with refs.||14|| | |25| | ||17|| . Then by ( |4.5| ), U(<p) is normalised as 



U"(0) = 1 . (4.9) 

Again, the conclusions of this section are not particular to the LPA, and can be expected 
to hold also for the exact RG, with one proviso: beyond the LPA we would have to consider 
the effects also of wavefunction renormalisation. This is straightforwardly incorporated by 
a further change of variables (multiplicative renormalization) to a field with a kinetic term 
that remains correctly normalised under flow in t (or to). As a byproduct one obtains a 
non-perturbative expression for the running anomalous scaling 7(771). We will provide the 
details in sect. 7, where the 0(d 2 ) approximation will be developed. 



5. Estimates from the LPA of the Wegner-Houghton RG 

In this section we carry out explicitly the program outlined above, for equation (|2.2|) , the 
result of the LPA applied to the RG for sharp cutoff. We will also explain our numerical 
methods for solving the resulting partial differential equations to obtain the universal 
coupling constant ratios i.e. the Taylor expansion coefficients of U((p). Since these results 
refer to the case of the Wilson-Fisher fixed point in three dimensions, we from now on set 
D = 3, implying for the LPA, d = 1/2. 

Our first step is to define an appropriate running mass to. If we set 

a(t) = V"(0,t) and (for later) &(t) = F (4) (0, t)/24 , (5.1) 

then by Q and flOD , 

a/m 2 -> 1 , (5.2) 
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as both tend to infinity. We cannot simply have m 2 = a however, since a starts out 
negative: 

a(t = -oo) =a* = -.46153372- •• (5.3) 



||22j| PI , this being tr's fixed point value. Using this in eqn.( [2.2| ), the Taylor expansion of 
V((f>) can be developed and hence of course fixed point values for a and all the other 



couplings may be deduced. Some of these were quoted in ref. ||28|| . Since there has been 
confusion in some of the recent literature, it is worthwhile to stress that these numbers 
have nothing directly to do with the universal ratios in the scaling equation of state. In 
particular, fixed point values depend sensitively on the form of the cutoff; they are not 
universal. 

Expanding (|2.2|) in powers of the field, and using (|5.1|), we have 



d° ~ 24a 

aF = 2ff + TT7 ' (5 ' 4) 



Clearly in general, on expanding (|2.2| ), inverse powers of 1 + a are generated, so to keep 
the algebra clean it is helpful to define 



1+m 2 , (5.5) 



in agreement with ( |5.2| ) . At the same time this ensures that m is real, since from ( |2.2| ) we 
clearly have a(t) > —1 for all t. Now define the couplings in U((p, m) by, 

oo 

U(ip,m) = ~—~ip 2 + U int (<p,m) with U int (<p, m) = £ + ^^(m) . (5.6) 

k=2 



Thus, from (|5.1j ) and (|4.7|) we have = a/m. Differentiating ( |4.7|) with respect to t, and 



using this relation, (|5.4|) , (|5.5|) and (|2.2|) , we obtain 



P^U(^m) = (^- 1 ^)Uu- 1 -^U')+^Hl + Ul' nt ) , (5-7) 
am \m z m J J \ 2 / m- 3 

with the beta function, 

f3(m) = m — 1/m + 12a ^(m) / m 2 . (5-8) 

The ip 2 part of ( |5.7| ) is satisfied automatically by (|5.5| ) , while all other couplings, including 
«4, are deduced as a function of m from ( p77|) once the boundary condition at m = m* = 
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\/l + cr* is provided. This boundary condition may be deduced as follows. First we need 
to choose a normalization for v(4>) (cf. sect. 3). We choose, 



w"(0) = l , (5.9) 

which by ( |4.1| ) implies for small g(t) (i.e. t — > — oo), a = a* + g(t), and thus at the fixed 
point dg(t)/dm = 2m*. Defining, analagously to (|4.7| ), 

U*((p) = V((p^/m^)/ml and u(ip) = v(cpy^mZ)/ml , (5.10) 



substituting (|4.1| ) into ( |4.7|) and differentiating with respect to to gives the boundary 
condition for dU/dm, 



dU 
dm 



— {±<pUl-3U*)+2m*u . (5.11) 
to* 



Note that dU /dm needs to be specified as well as U at the fixed point, because /3(to*) = 
and the right hand side of (|5.7|) vanish there. [Thus m = to*, 11 = 11* corresponds to a 
so-called 'singular point' for the differential equation (|5.7|)-1 

Since our goal is the numerical computation of as many of the universal coupling 
constant ratios a2A;(oo) [k = 2,3, • • •, cf. eqn. fl4.5|) and below (|4.7|) 1 as possible, it makes 
sense to Taylor expand the differential equation (|5.7|) and boundary condition (|5.11|) in 
terms of (p and obtain directly ordinary differential flow equations for the a,2k{iri)- These 
are easily obtained to high order using algebraic computing packages. (We used Maple.) 
We can expect to obtain better estimates this way and with less effort compared to, say, 
numerically solving the partial differential equations for U((p, to), and then obtaining the 
ct2fc(oo) by numerically differentiating U(<p,oa), 2k times at (p = 0. But we now have 
a problem of truncation: the expansion of ([5.71) leads to differential equations where the 
da2k/dm depend on all «2j up to and including «2fc+2, and therefore if we keep only a 
finite number of these equations, say up to and including the equation da2k/dm = • • •, we 
will require an ansatz for the highest <^>-power, the 'top coupling' «2fc+2(^)- However, this 
problem is not as severe as for computations of fixed points by truncations in powers of 
the field||, in particular since we already know the initial (fixed point) data numerically 
there will be no spurious solutions @, and we can expect that the results converge as we 
keep more and more of the equations (and thus at the same time compute estimates for 
more and more of the couplings) providing we supply a decent ansatz for the top coupling. 
Of course if our ansatz corresponded exactly to the right answer for the top coupling, 
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all the lower couplings would be computed correctly. Our assumption is however that if 
a reasonable model of the top coupling's flow is inserted, then the lower couplings flow 
becomes more and more insensitive to the error in the modelling of the top coupling as k 
is increased (the error having to feed down successively through more equations) . We find 
that the results bear out this assumption for several different ansatze, at least to sufficient 
accuracy for this exercise. These ansatze are outlined below. First we sketch the final 
steps necessary to obtain appropriate boundary conditions for these flow equations. 

From the numerical value for cr*, we deduce as many of the couplings in V(<fi) as needed 
(c.f. ( |5.3| ) and below it), and thus from ( |5.5| ) and ( |5.10| ), as many of the initial values, 
«2fc(w*), as needed. Similarly, by Taylor expansion of (|3.2|) , using the normalization ( |5.9|) , 
and the relevant eigenvalue A = 1.450416(1) (already obtained numerically in the process 
of writing refs.@@) we can obtain the Taylor series coefficients of v(<p). From ( |5.10| ) 
and ( |5.11|) , this then gives the numerical values for the initial gradients a' 2k (rn*). Some 
numerical methods can cope with this singular point [c.f. below (|5.11| )], but we found that 
the resulting differential equations were remarkably well-behaved numerically and that 
such sophistication was not necessary. Instead we used the basic fourth-fifth order Runge- 
Kutta provided within the Maple package, and set the boundary condition a little way 
away from the fixed point as a^A^ra* + S) = a.2k(jn*) + Q! 2fc (m*) Here, 5 could be set 
very small. We chose typically to set it to ~ 10 -6 . The results were entirely insensitive 
to the precise choice. We integrated out to m = 20: with the improved estimate ( |5.13j ) 
described below this gave us more than sufficient accuracy while being easily achieved by 
the above package. 

Turning now to the ansatz for the top coupling, the simplest workable ansatz of all, is 
to ignore its evolution and just set it equal to its initial value: 



Numerically we found that this ansatz results in reasonable convergence, but it is easy 
to do better. From the above, we know also the initial gradient, while from the previous 
analysis (and borne out e.g. by the above ansatz) we expect that all couplings a,2j(m) 'run' 
with m initially and then freeze out at some value mf reeze , after which a2j(m) ~ a!2j(oo). 




(5.12) 



Indeed from ( |5.7| ) we see that asymptotically ot2j(m) ~ (1 
used to improve the estimate of 0^ (00): 



3-fc 
2m 2 



)a2j(oc), a fact which we 




(5.13) 
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This suggests the following simple ansatz for the top coupling, 

a 2 fc+2(m) = a 2 k+2{m Sf ) H r r at 2k+2 (m*) , (5.14) 

which thus has the right value and gradient at m = m* and will freeze out at the same 
value mf reeze as a^. Initially we employed an ansatz where «4 was replaced in (|5.14[) by 
«2fc, in the expectation that any k dependence of m/ reeze would thus be better taken into 
account: 

/ \ / \ i at 2 k(m) - a 2k (m*) , 

«2fc( m *) " " 

It is such an ansatz that led to the 0(d 2 ) estimates quoted in ref. [14||. However, we found 



that there was little or no drift of m/ ree2e with k, while using a 2 k in place of «4 in ( |5.14| ) 
led sometimes to numerical instabilities e.g. when ct' 2k (m^) happened to be much smaller 
than «2fc+2( m *)' an d to poorer convergence. 

All three methods described led to consistent estimates for the a2fc(oo). We present 
some of the raw data in sec. 8. The most convergent answers were produced with ansatz 
( |5.14| ). Using these, and untying the numerical transformation ( |2.4| ), i.e. as 

«2 fc = C 2fc - 2 « 2fc (oo) , (5.16) 

where by below ( |2.4| ), £ = 2n, yields our final results. These are given in table 7. They 
correspond to the normalization used in ref. (TTj . 



6. Lowest order in the derivative expansion 

In this section we carry out the program for the LPA of the RG written directly in 
terms of the Legendre effective action, c.f. (|2.5|), with respect to a smooth cutoff. This 
corresponds to O(d ) in a derivative expansion of T^: 

r A [0] = Jd\ {i(a^ phys ) 2 + y phys ((/) phys ,A)} , (6.i) 

with no other approximation made. (We drop the c superscript.) The relevant flow equa- 
tions were derived in ref. 0, where details of the form of the cutoff etc. can be found. 
Since the critical exponents were found [0 to be better estimated by the smooth cutoff 
LPA compared to the sharp case described in the previous section, one might expect the 
universal couplings also to be better estimated. The results bear out these expectations. 
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Instead of the flow equation (j2.2|) , one finds 0: 

d_ 
dt 



—V((f>, t) + \<j>V'{<j>, t) - 3V(</>, t) = -l/y/2 + V"(c/>, t) . (6.2) 



We define again a(t) and a(t) by eqn.(|5.1j). However, Taylor expanding ( |6.2| ) one finds 
instead of ( |5.4|) , 

^ = 2a + 12 & (2 + a)" 3 / 2 , 

the difference of course arising as a result of the inherently different 'blocking' procedure 
implied by use of the sharp cutoff. (Only universal quantities such as the exponents and 
the a2k(oo) are independent of the choice of cutoff and this only when the RG is computed 
exactly, otherwise some dependence on the choice of cutoff is left||. See also our comments 
below Q5.3Q .) By similar reasoning to that for (|5.5|) , a convenient definition for m is here, 



a = -2 + m 2 . (6.3) 
Using the same expansion (|5.6| ), we thus obtain the flow equation 

with the beta-function, 

f3(m) = m — 2/m + Qa^{m)/m 2 ' 



Choosing the normalization (|5.9| ) for the relevant perturbation, and using Qf.l|) , 
again implies dg(t)/dm = 2m* as t —>■ — oo, and thus with definitions (|5.10| ), the same 
formula ( |5.11| ) holds for the initial boundary conditions. Our treatment of these and the 
choice of ansatz for the top coupling are the same as in the previous section. From ( |6.4| ) , 
the analogous formula to ( |5.f 3| ) reads however, 

ajyCoo)- (l + ^ja 2i (m) + oQjJ . (6.5) 

The numerical values for the initial couplings follow from[[| a* = —.534648257 • • •, and the 
relevant eigenvalue A = 1.514260- • •. Eqn. (|6.2| ) was also prettified by a numerical trans- 
formation of the form (|2.4j), and thus must be untied as in (|5.16| ), but the corresponding 
value for ( is £ = V2tt 0. Numerical results for a^ k are given in sec. 8. 
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7. Order derivative-squared 

We use the smooth cutoff Legendre flow equation of the previous section, but improve 
the approximation by keeping aiso ail terms of 0(<9 2 ), i.e. we replace (|6.1|) with 



f a[</>] = Jd 3 x {^(d^ phys ) 2 K phys ((j) phySl A) + "l/ p hys(</>phys,A)} ■ (7.1) 

No further approximation is made. The resulting flow equations were derived in ref. || : 
dV 1 , l-r?/4 

™ + i(l + r)W -3V= - % = (7.2a) 

dt 2 \[K y/V" + 2VK 



and — + -(1 + riUK' + riK = (l-~) ^ ; >- ; - 7.26 

dt 2 K ' V 4/ [ 48 K* 2 {V" + 2VK) s 2 



1 58V"'K'VK + 57(K') 2 + (V"') 2 K 5 (y'") 2 ^ + 2V"'K'VK + (K 



l\2 



+ 



48 K(V" + 2VK) 5 / 2 12 v / i?(^" + 2v / if) 7 / 2 

The main novelty in the theoretical construction, compared to the previous sections, is 
the appearance of a non-zero anomalous dimension 77. This is determined at the fixed 
point itself, as consequence of a rescaling invariance which is here preserved by the 
derivative expansion as a result of the careful choice of cutoff function. We refer the reader 
to refs.||[|| for the details. They will not be important for the ensuing analysis. The 
non-zero value for 77 means that the couplings have anomalous scaling dimensions. These 
can be taken into account by approprate insertions of m v , however, with an eye to future 
more general applications, we will choose the more conventional approach and introduce 
wavefunction renormalisation. 

The change to dimensionless quantities is given, c.f. (|2.3|) and ref. [EL as 



<P = <P P ^ 2 /A d , V(<p,t) = V p (<P pi A)/A 3 and K(cj>,t) = (A/^K^A.) , (7.3) 

where now d = |(1 + 77), and here the powers of \i v are introduced to balance engineering 
dimensions. These physical quantities are related to the true physical quantities by the 
numerical transformation fl2.4p , and K p = K p h ys p|, with ( = \/2~n as in sec.||. As 
previously, we will correct for this transformation at the end of the section. 

The physical fields do not yet have a normalised kinetic term. It is helpful to define 



z(t) = l/VK(0,t) 
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(7.4) 



then from (|7.3j ) and ( |7.1| ), the renormalised physical fields 



(f>R = e^^/zit) (7.5) 

have a normalised kinetic term Kb.{4>r-, A) = K p ((f> p , A)/K p (0, A), satisfying Kr(0, A) = 1. 
(The normalisation of Vr(4>r, A) = V p ((f)p, A) is of course unaffected.) Now we can define 
the physical mass scale [c.f. (fOf)], 



mJ = V&'(0,0) , (7.6) 

and thus the universal dimensionless potential and 0(p 2 ) parts, 

U(<p) = V R (<p y /m;,0)/m 3 p , L{<p) = K R (<py/m;,0) ■ ( 7 - 7 ) 

Note that U(f>) is still normalised as in Q4.9|) , while evidently L is normalised as 

L(0) = 1 . (7.8) 

Introducing again a running mass m(t) with the property that it satisfies the limit ( f4.6|) , 
and combining the above transformations, we see that the quantities 

U((p,m) = V((pzy/m : t)/m s and L((p,m) = z 2 K(ipzy/m,t) , (7-9) 

thoroughly enjoy the property that they tend to their universal counterparts as m — > oo. 
Defining again a{t) as in ( |5.1|) , we see from ( |779|) and (^4.9|) , that (|5.2|) is replaced by the 
condition 

az 2 /m 2 -> 1 . (7.10) 

Since K((f>,t) starts out normalised at the fixed point, i.e. K(0, — oo) = K(0) = 1 0, we 
have by (|7.4[) , the initial condition z(t = — oo) = 2;* = 1. On the other hand a(t = —00) = 
a* and we found numerically (in the process of writing ref. [H), a* = —.3781684 • • -, so once 
again replacing the limit by equality in (|7.10| ) will not do. By a similar analysis to that 
leading to ( |4.3| ) and (O), we may confirm that K p ((j) p ,A) has a finite limit as A — ► 0, 



and thus determine from (|7.3j ) that z ~ m v ^ 2 as m — >• 00. Since 77 is numerically small, 
77 = .05393208 • ■ ■ @, the choice 

m 2 =az 2 + 2z (7.11) 
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satisfies ( |7.10| ) while, being proportional to a + 2/z, neatening up the small field expansion 
of eqns(|7.2|). Now, defining the couplings similarly to (|5l3), 



1 az 2 °° 
U(<p, m) = - — T ip 2 + U int (<p, m) .Uintiip.m) = £ + y^ip 2k a 2 k(m) , (7.12a) 



2 m 2 

k=2 



L(<p,m) = 1 + J> 2fc 7 2fc(™) , (7.126) 



k=i 



and substituting (|7.9| ) into eqns( |7.2[ ), we obtain 



P-£-Ufa m) = ^( \pU> - 3U) + - (l - — (7.13a) 

„ 9 Tl v , x A r/ r \ r/ / r]\ ( z 2 24LL" -19(L') 2 . 

~dm ™> = ^ " (2^ + £ ) + ^ £ + (' " j) {48 n,W/V (™») 

2 58m 2 zU'"L'VL + 57z 2 (L') 2 + m 4 (U'") 2 L + bz 2 m 4 (U"') 2 L + 2m 2 zU"'L'VL + z 2 (L') 2 



48 m 6 i? 5 L 12 m 8 R 7 VZ 

where we have introduced the short-hand, 



and defined the beta function /3(m) = dm/dt and running anomalous dimension, 

7(m) = iat • 

The factor of 2 here corresponds to the usual definition in terms of a replacement z = \f~Z 



in (|7.4f) and ( f7.5|) . These functions may be derived analagously to that of (|5.4| )-(p7q), or 
directly by evaluating the tp 2 coefficients in ( |7.13o| ), and the constant term in (|7.13£| ), using 
( gTg ). Thus, 

V2 / m V 4/ \ m 13 m° / ,^ ^ 

To determine the initial boundary conditions we again use ( [4.1| ), together with the K 
component: 

K(<l>,t) = K(<f>)+g(t)k(<f>) , (7.15) 
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where k(<p) is the 0{d 2 ) part of the relevant operator 0. It is simpler to implement the 
normalisation k(0) = 1, which thus replaces ( |5.9|) . If we define for convenience, r = v "(0), 
then (U), ( fTT5| ), ( PD and d774|) imply that for g(t) small, i.e. t ^ -oo, 



a(t) = (7* + Tfif(t) 



and 



l/z 2 (t) = l + g(t) 



Substituting (|7.11|) and differentiating with respect to m we deduce that at the fixed point, 

dg _ ^dz _ 2m* 
dm dm 1 + r — 

Thus, using the same definition ( |5.10[ ) (because z* = 1), and analogously 

L*(ip) = K{tp^[ml) , l((f) = k(ipy^m^) , 



combining ( |4.1[ ), ( |7.15[ ), and (|7.9|), and differentiating with respect to m gives the boundary 
conditions for the gradients, 



dU 




1 


dm 


* 


m* 


dL 




1 


dm 


* 


2m* 



2 



m; 



1 -r 



m* — 1 — r 
2L + <pL' - 21) 



(<pU' - 2u) 



As before, we perform a small-field expansion of these equations, determining the initial 
data for the couplings in ( |7.12| ), in the way described in sec.[5], using the numerical values 
of cr* and r\ already given, and (from the analysis that led to ref.0) A = 1.61796660 ■ ■ •, 
r = -24.472671 • • •. 

Expanding ( |7.13|) in powers of cp, one sees that the flow equations for the couplings 
have the following structure: two infinite sequences of equations, the differential equations 
for «2fc, with k = 2, 3, • • •, and the 72/c equations with k = 1, 2, • • •, each of form, 

d 

— a 2 k(m) = / 2 a fc(«4, ■ • • , a 2 k+2] 72, • ■ • , 72fc) , 



_d_ 

dm 



l2k{m) = fl k {&± 



. 72fc+2) 



(for some functions f^., f^k)- This structure implies that, for given k > 2, there is not one 
unique minimal way to close the equations, rather there are two. These are: 
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i) The differential equations for c*4, • • • , a 2 k and 72, • • • , 72fc-2 are kept, with ansatze being 
supplied for 7 2fc and a 2k+2 . 

ii) The differential equations for a^, • • • , a 2k and 72, • • • , 72fc are kept, with ansatze being 
supplied for 72^+2 and a 2fc+ 2- 

Since there is no a priori reason to prefer one method over the other, we use both. 
Indeed, by employing at each k, first method (i) and then method (ii), we increment 
sequentially by one each time, the number of couplings which are estimated. For those 
top couplings that require an ansatz we continue to use ( |5.14j ), with also a 2k + 2 replaced 
by 72/c or 72/C+2 as appropriate. (Of course other ansatze are possible.) 

From a large m analysis of ( [7.13|) and (|7.14| ) , we determine that an improved estimate 
for «2j(oo) can be obtained as 

(z (ffl) \ 
1 + (3- j)—T J a 2j(m) , 

which thus replaces ( |6.5| ) . The 72 j have also z / m 2 corrections but with a more complicated 
coefficient [the result of expansion of (U"') 2 (l + f7") _5/2 ]. Rather than compute this 
coefficient algebraically, we determined it by a simple numerical fit on neighbouring values 
of m in the solution 72j(m). This gave more than sufficient accuracy for the improved 
estimate. 

Finally, we have to fold back in the missing numerical factors [c.f. below ( |7.3| )1, thus 
the final results are given as 

a p 2k = ( 2k - 2 a 2k (oo) and 7^ = ^72.(00) , (7.16) 

where ( = \^2n. 



8. The numerical results 

In this section we display the numerical results obtained from the analysis described in 
the previous three sections. In table 1, we compare the results obtained for a^, at O(d ) 
approximation as in sec.[| employing the three different closure ansatze, ( |5.12| ) , ( |5.14| ) and 
( |5.15| ) for the 'top' coupling a 2 k+ 2 (m). One sees that indeed there is apparent convergence 
for each method as k is increased. We provide also in table 1 an estimate of the converged 
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answer for each ansatz, together with an error which reflects the spread of results over the 
several largest values of k in the table. 



2 
3 
4 
5 
6 
7 
8 
9 

10 



(EUD 

1.1955 
1.2917 
1.3117 
1.3066 
1.3004 
1.2994 
1.3008 
1.3017 
1.3016 



(EH) 
1.3137 
1.3097 
1.3052 
1.2954 
1.3003 

1.3017 
1.3014 



(gg) 

1.3137 
1.3097 
1.3041 
1.3008 
1.3003 
1.3010 
1.3014 
1.3014 
1.3012 



oo 



1.3016(8) 



1.301(1) 



1.3012(2) 



Table 1. O(d ) results for a^, at different levels k of truncation, using three different 
closure ansatze. For the two missing values in the middle column, the integrating routine 
failed to converge, due to wild numerical behaviour of the ansatz [c.f. below (|5.15[) -1 

The three methods give consistent estimates, although method ( |5.14| ) is seen to be the 
most powerful. Since we find that these facts hold true for the other couplings (and the 
other flow equations), from now on we present numerical results for method ( |5.14| ) only. 



k 


a\ 


«6 




"lO 


"l2 


"14 


"16 


"18 


«20 


2 


1.3137 


















3 


1.3097 


2.7808 
















4 


1.3041 


2.6873 


1.964 














5 


1.3008 


2.6556 


1.704 


-3.566 












6 


1.3003 


2.6537 


1.730 


-2.718 


-1.83 










7 


1.3010 


2.6600 


1.779 


-2.292 


2.05 


17.1 








8 


1.3014 


2.6632 


1.795 


-2.231 


2.00 


9.72 


20.0 






9 


1.3014 


2.6628 


1.790 


-2.286 


1.46 


4.70 


-24.8 


-47.3 




10 


1.3012 


2.6615 


1.783 


-2.324 


1.25 


3.73 


-26.4 


-16.9 


-389 


oo 


1.3012(2) 


2.662(2) 


1.78(1) 


-2.3(1) 


1.3(7) 


6(4) 


-20? 


-20? 


? 



Table 2. O(d ) results for all computed couplings ckL, at different levels k of truncation 
via ( |5.14| ), together with estimates for k = oo. 
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k 




a? 

D 


at 

o 


1U 


a? 9 

1Z 


2 (i) 


.8682 










2 (ii) 


.8688 










3 (i) 


.8644 


1.253 








3 (ii) 


.8667 


1.269 








4 (i) 


.8639 


1.237 


.586 






4 (ii) 


.8646 


1.240 


.5790 






5 (i) 


.8639 


1.235 


.5587 


-0.950 




5 (ii) 


.8633 


1.230 


.5275 


-1.112 




6 (i) 


.8635 


1.232 


.5475 


-0.873 


-.176 


6 (ii) 


.8630 


1.229 


.5359 


-0.899 


-.0688 


oo 


.8635(5) 


1.230(3) 


.54(1) 


-.95(5) 


? 



Table 3. 0(d 2 ) estimates of a\ - at different levels of truncation k via (|5.14|) , for the two 
closures described at the end of sec. [7], together with estimates for k = oo. 



k 


7 2 P 


74 P 


7e P 


7 8 P 


7 P o 


7l2 


2 (i) 


.6427 












2 (ii) 


.6472 


-3.023 










3 (i) 


.6442 


-3.047 










3 (ii) 


.6537 


-2.748 


7.26 








4 (i) 


.6504 


-2.782 


7.10 








4 (ii) 


.6524 


-2.771 


6.59 


7.82 






5 (i) 


.6514 


-2.776 


6.55 


7.3 






5 (ii) 


.6499 


-2.801 


6.05 


-2.5 


-69 




6 (i) 


.6502 


-2.797 


6.10 


-2.20 


-71.3 




6 (ii) 


.6492 


-2.805 


6.02 


-2.4 


-53.7 


-15 


oo 


.6497(5) 


-2.805(5) 


6.06(4) 


-2.3(1) 


-62(9) 


? 



Table 4. 0(d 2 ) estimates of 7^ at different levels of truncation k via ( |5.14j ), for the two 
closures described at the end of sec. [7], together with estimates for k = 00. 

The estimates resulting from LPA of the Wegner- Houghton equation are given in table 
5. (The accuracy of these numbers was improved by correcting also for the 1/m 3 contri- 
butions in ( |5.13| ), in a similar way to that done for 7^ in the previous section.) To our 
knowledge these estimates for universal coupling constant ratios have not been derived 
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before. (N.B. this approximation gives v = .6895 and uj = .5952 f2j| ||22|| .) 



k 


a? 

4 


a? 




a? 


1U 


a? 9 

LZ 


14 


10 


1© 


2 


1.5308 
















3 


1.5288 


3.569 














4 


1.5184 


3.439 


2.142 












5 


1.5107 


3.377 


1.833 


-5.180 










6 


1.5122 


3.396 


2.005 


-3.491 


5.88 








7 


1.5163 


3.426 


2.142 


-2.944 


7.30 


14.4 






8 


1.5158 


3.421 


2.094 


-3.317 


4.48 


-5.50 


-58.6 




9 


1.5135 


3.404 


2.027 


-3.556 


3.63 


-6.49 


-36.3 


31 


oo 


1.514(2) 


3.41(1) 


2.08(6) 


-3.3(3) 


5(2) 


-6? 


-40? 


? 



Table 5. Sharp cutoff (Wegner-Houghton) with LPA, at different levels k of truncation via 
( |5.14|) , together with estimates for k = oo. 

We can attempt to combine the O(d ) and 0(d 2 ) results of tables 2 and 3, to give 
numbers with an error which takes into account that due to truncation of the derivative 
expansion. Although two terms in the derivative expansion approximation are not really 
enough to take this exercise too seriously, we note that for all the results, including expo- 
nents 11,0 the accuracy of the numbers obtained (compared to results from more accurate 
methods where these are available e.g. in ref. ||14|| ) , obey the following pattern: the sharp 
cutoff results are the worst, O(d ) are better, and 0(d 2 ) better still. Comparing the O(d ) 
and 0(d 2 ) estimates of exponents^ v and u to the worlds best determinations we find 
that 1/3 of the difference between the O(d ) and 0(d 2 ) results give a good estimate of the 
error in the 0(d 2 ) results. If we adopt this algorithm for the couplings a^- we obtain the 
numbers given in the next table. We include the data from ref. , and the most accurate 
estimates for the a^- - three dimensional resummed perturbation theory - reported by 
ref. |TJ|] [translated to this normalisation through eqn.(8.2)]. 

We have not included in table 6, the 0(d 2 ) estimate of rj [see above ( [7.11 )1, or the 
0(d ) estimates of a\ 2 and given in table 2, because we have only one term in the 
derivative expansion for these and so cannot estimate the error due to truncation of the 



14 We stress again that these do not correspond to exact RG fixed point couplings [c.f. our 
discussion above (5.4)]. 



15 For the sharp cutoff case see above, and refs. |2[ p^] 

16 The 0(d 2 ) estimate for uj in table 6, corrects a misprint in ref.Q. 
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Approx' 


V 






a 6 


a 8 


a 10 


o(d°) 

0(d 2 ) 


.6604 
.6181 


.6285 
.8972 


1.3012(2) 
.8635(5) 


2.662(2) 
1.230(3) 


1.78(1) 
.54(1) 


-2.3(1) 
-.95(5) 


d exp n 


.618(14) 


.898(90) 


.86(15) 


1.2(5) 


.5(4) 


-1.0(5) 


mi & u\ 


.631(2) 


.80(4) 


.988(2) 


1.60(1) 


.83(8) 


-2.0(1.3) 



Table 6. Estimates from the derivative expansion, with errors computed as described in the 
text, compared to combined results from the worlds best determinations for the exponents 
29||P1, and resummed perturbation theory for the a p 's [fT4| |. The line labelled "<9 exp n " 



gives a combined result from 0(d ) and 0(<9 2 ), as explained in the text above. 



derivative expansion. Similarly, at 0(d 2 ) we obtain for the first time, estimates of 7^, 
as displayed in table 4. Note that the truncation level where we choose to stop in each 
of these tables, was determined by limitations of size and/or stability within the Maple 
computing package. It is certainly possible to do better, and thus estimate more higher 
order couplings, with a more serious attack e.g. by writing FORTRAN code, nor do we 
wish to rule out the possibility that a method based on direct integration of the partial 
flow equations derived in the previous sections, could produce competitive or better results 
than obtained here (while having the advantage of course of obtaining directly also the 
equation of state). 

It can be seen from table 6, that while perturbative methods are more powerful than 
the derivative expansion for low order couplings, the derivative expansion eventually wins 
out. The reason for this is that the derivative expansion at these lowest orders, is crude in 
comparison to the perturbation theory methods, however the perturbative methods suffer 
from being asymptotic - which in particular results in rapidly worse determinations for 
higher order couplings. The derivative expansion does not suffer from this, since it is not 
related at all to an expansion in powers of the field. Indeed, it may be shown that|| even at 
the level of the LPA, Feynman diagrams of all topologies are included. We recall here that 
derivative expansion estimates can successfully be given for the multicritical fixed points 
in two dimensions ||, where all other standard methods fail.0 In that study, we found that 
derivative expansion estimates actually tend to improve for higher dimension corrections 

Needless to say, it is possible, and would be very interesting, to apply the methods developed 
here to the two-dimensional cases and thus derive scaling equations of state for each of these 
multicritical points. 
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to scaling (i.e. higher dimension operators) and for greater multicriticality. This trend 
is opposite to that in perturbation theory (or any other standard approximation method 
for that matter). This suggests that the main source of error in the derivative expansion 



estimates of the a^/s is that of the coupling. We follow Guida and Zinn- Justin [|~4 
and factor out the a\ couplingjl! by defining 

f(z) = 24a 4 (oo) [U (z/ v / 24a 4 (oo)) - 

= z 2 /2 + Z y4\ + J2z 2l F2i-i/2l , (8 - 1} 

1=3 

c.f. OOP , (U) and (|3). The F's are then given by 

F 2 ;_ 1 = 2/a 2Z (oo)/[24a 4 (oo)]'- 1 = 21 a p 2l /(24 al) l ~ l . (8.2) 

Note that they are independent of the numerical normalisation factors (|7.16| ) . We find much 
improved derivative expansion estimates of the higher order F's compared to the ct p 's. This 
is illustrated in table 7. It seems therefore that the main source of error in the derivative 



expansion can be absorbed in an effective normalisation factor ( e ff = Q^J (a^) exact /a^. At 
O(d ), Ceff = 1-148(1) C At 0(d 2 ), Ceff = .935(1) C- We have taken the opportunity to 
include in table 7, the estimates from the sharp cutoff LPA data in table 5, and comparisons 
with estimates from other methods.0 We use the same heuristic as before to combine the 
0(<9 ) and 0(d 2 ) results into a single derivative expansion estimate with error. In addition, 
at O(d ) we obtain the estimates F\\ = 5(3) x 10~ 7 , F i3 = 9(6) x 10~ 8 , and with sharp 
cutoff LPA, the estimate Fn = 9.5(4.0) x 10~ 7 (consistent with the O(d ) result). To our 
knowledge there are no other estimates with which these can be compared. 



18 Ref. [|4}'s coupling g* = 2Aa\. 

19 We recall again that the 0(<9 2 ) estimates quoted in the corresponding table in ref. [14] were 



derived using the less powerful closure ansatz ( |5.15 ). They are consistent with the present es- 
timates but with a larger error - reflecting the poorer convergence with truncation level k. For 



ref.[17], we list only their more sophisticated Pade-Borel estimates. 
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Approx' 
Sharp 

o(d°) 

0{d 2 ) 



1.514(2) 

1.3012(2) 

.8635(5) 



.0155(3) 

.01638(1) 

.01719(4) 



3.6(5) x 10" 4 
4.68(3) x 10" 4 
4.9(1) x 10" 4 



^9 

-1.7(5) x 10" 5 
-2.4(1) x 10" 5 
-5.2(3) x 10" 5 



d exp n 



.86(15) 



.0172(3) 



4.9(1) x 10" 4 



-5(1) x 10" 



4.96(49) x 10" 4 
4.1 - 16.2 xlO" 4 
4.5(3) x 10" 4 

4.3 x 10" 4 



5.4(6) x 10" 4 
.00236(40) 



d = 3 
d = 3 

e— exp. 

e— exp. 





m 

erg m 

HT [gO | 
HT [[!T| 
HT g2| 
MC p| 
MC p| 



.988(2) 

1.2 

1.20 

.99(6) 

1.019(6) 

.987(4) 

.97(2) 

1.020(8) 



.01712(6) 
.0168 - .0173 
.0176(4) 
.0176 
.016 

.0205(52) 

.01780(15) 

.017(1) 

.0227(26) 

.027(2) 



-6(4) x 10" 



-3.2(2) x 10" 



-2(1) x 10" 



Table 7. Estimates for a\ and the i*2Z-i : from the derivative expansion - with errors 
computed as described in the text, the LPA of the Wegner- Houghton equation, and (sum- 
maries of) other methods. Again the line labelled "9 exp n " gives combined results from 
0(d ) and 0(9 2 ), in a way described earlier in the text. 
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